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ABSTRACT 

We correct the observed Milky Way satellite luminosity function for luminosity bias using published com- 
pleteness limits for the Sloan Digital Sky Survey DR5. Assuming that the spatial distribution of Milky Way 
satellites tracks the subhalos found in the Via Lactea ACDM N-body simulation, we show that there should 
be between ~ 300 and ~ 600 satellites within 400 kpc of the Sun that are brighter than the faintest known 
dwarf galaxies, and that there may be as many as ~ 1000, depending on assumptions. By taking into account 
completeness limits, we show that the radial distribution of known Milky Way dwarfs is consistent with our 
assumption that the full satellite population tracks that of subhalos. These results alleviate the primary wor- 
ries associated with the so-called missing satellites problem in CDM. We show that future, deep wide-field 
surveys such as SkyMapper, the Dark Energy Survey (DES), PanSTARRS, and the Large Synoptic Survey 
Telescope (LSST) will deliver a complete census of ultra-faint dwarf satellites out to the Milky Way virial 
radius, offer new limits on the free-streaming scale of dark matter, and provide unprecedented constraints on 
the low-luminosity threshold of galaxy formation. 

Subject headings: cosmology: observation — surveys — galaxies: halos — galaxies:dwarfs — galaxies: Local 
Group 



1. INTRODUCTION 

It is well-established from simulations in the ACDM con- 
cordance cosmology that gala xy halos are forme d by the 
merging of smaller halos (e.g. IStewart et al.l 120071 and ref- 
erences therein). As a result of their high formation red- 
shifts, and correspondingly high densities, a large number 
of self-bound dark matter subhalos should survive the merg- 
ing process and exist withi n the dark matter halos of 
galax i es like the Milky Way dKlypin et alJl 19991: iMoore et aT] 
119991: IZentner & Bullockl |2003|) . A direct confirmation of 
this prediction h as yet to occur. Surveys of the halos of the 
Milky Way (e.g.lWil lman et a f]|2005t IBelokurov et alj|2006t 
IWalsh et all l2007at IBelokurov et al.N2008j) and Andromeda 
(e.g. iMartin et all 120061: iMaiewski et all 1200H llrwin et all 
l2008tlMcConnachie et aT]|2008l) have revealed only ~ 20 lu- 
minous dwarf satellites around each galaxy, approximately 
1 order of magnitude fewer than the expected number of 
subhal os that are thought to be massive enough to form 
stars dKlypin et al.l 11999c IMoore et al.l [l999t iDiemand et all 
l2007bi IStrigari et al.l l2007bl) . Of course, the mismatch 
is much more extreme (larger than a factor of ~ 10 10 ) 
when compared to the full mass function of CDM sub- 
halos, which is expected to rise as ~ 1/M down to the 
small-scale cluste ring cutoff for the dark matter particle , 
M„,„ <g lMrr. (lHofmannetal.ll200lUBertschinged 12001 
Profumo et al. 2006; Johnson & Kamionkowski 2008) 

Astrophysical solutions to this missing satellites problem 
(MSP) include a reduction i n the ability of small halos to ac- 
crete gas after reionizatio n dBullock et alJ l2000: So merviliel 
l2002t iBenson et aT]|2002l) . and tidal stripping, which shrinks 
the m ass of halos after they have formed stars dKravtsov et al.l 
2004). Although there are some promising new tech- 
niques in development for studying t he formation of low- 
mass galaxies within a CDM context ( Rober tson & Kravtsovl 
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2007; Kau fmann et alj2007tlOrban et alj2008h . current mod- 
els have trouble explaining many details, particularly re- 
garding the lowest luminosity dwarfs (see below). Other 
more exotic solutions rely on dark m a tter particles that 
are not cold dHogan & Dalcantonl l2000t iKaplinghal 120051: 
Cembranos et al. 2005; S trigari et alJl2007d) or non-standard 
models of inflation (IZentner & Bullocki2 003), which produce 
small-scale power-spectrum cutoffs at M cut ~ 10 6 — 10 8 M Q , 
well above the cutoff scale in standard ACDM. In principle, 
each of these scenarios leaves its mark on the properties of 
satellite galaxies, although multiple scenarios fit the current 
data. Hence, precisely determining the shape and normaliza- 
tion of the lowest end of the luminosity function is critical to 
understanding how faint galaxies form and how the efficiency 
of star formation is suppressed in the smallest dark matter ha- 
los. 

The Sloan Digi tal Sky Survey (SDSS; 

lAdelman-M cCarthv et all 120071) has revolutionized our 
view of the Milky Way and its environment, doubling the 
number of known dwarf spheroidal (dSph) galaxies over 
the last s everal years (e. g. Willm an et alj|2005c IZuckert 
2006b Tal: IBelokurov et al ]|2006c iGrillmairl |2006t TWalsh"! 
2007a; llrwin et alll2007l) . Many of these new dSph galaxies 
are ultra-faint, with luminosities as low as ~ 1000L Q , faint 
enough to evade detection in surveys with limits sufficient 
for detecting most p reviously known Local Group dwarfs 
(IWhiting et al.1 120071) . In addition to providing fainter 
detections, the homogeneous form of the SDSS allows for 
a much better understanding of the statistics of detection. 
Unfortunately, given the inherent faintness of the newly 
discovered dSph's and the magnitude-limited nature of SDSS 
for such objects, a derivation of the full luminosity function 
of satellites within the Milky Way halo must include a 
substantial cor r ection for more distant undetectable satellites. 
Kopos ovet all d2007l) provided an important step in this 
process by performing simulations in order to quantify the 
detection limits of the SDSS and estimated the luminosity 
function by applying these limits to some simple radial 
distribution functions, finding ~ 70 satellites and a satellite 
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luminosity function consistent with a single power law of the 
formdN/dLv ~ L~ 125 . 

Our aim is to take the detection limits for SPSS Da ta Re- 
lease 5 (DR5) as determined by Kop osov et al.l d2007t) and 
combine them with a CDM-motivated satellite distribution. 
In order to provide a theoretically-motivated estimate for 
the total number of Milky Way satellite galaxies, we adopt 
the distributions of subhalos in the Via Lactea simulation 
dDiemand et alj|2007al) . and assume these to be hosts of satel- 
lite galaxies. Note that for the purposes of this paper, we de- 
fine a "galaxy" as a stellar system that is embedded in a dark 
matter halo. 

The organization of this paper is as follows: In $2] we de- 
scribe the overall strategy and sources of data used for this 
correction, as well as discussing the validity of the assump- 
tions underlying these data. fj3] describes in detail how the 
correction is performed and presents the resulting luminosity 
functions for a number of possible scenarios, while 2] dis- 
cusses the prospects for detecting the as-yet unseen satellites 
that the correction predicts, [^discusses some of the caveats 
associated with the technique used for this correction, as well 
as the cosmological implications of the presence of this many 
satellite galaxies. In [J6]we draw some final conclusions. 

2. APPROACH & DATA SOURCES 

2.1. Strategy 

The luminosity bias within the SPSS survey can be approx- 
imated by a characteristic heliocentric completeness radius, 
i? comp (Mv), beyond which a dwarf galaxy of absolute mag- 
nitude My cannot be observed. This radial incompleteness 
is accompanied by a more obvious angular incompleteness - 
the SPSS PR5 covers only a fraction / DR5 = 0.194 of the 
sky, or a solid angle JIdrs = 8000 square degrees. For the 
corrections presented in §3, we assume that magnitude My 
satellites with helio-centric distances r > R CO mp(My) have 
not been observed, while satellites with r < i? C omp(^v) 
have been observed, provided that are situated within the 
area of the sky covered by the SPSS footprint. Given an 
observed number of satellites brighter than My within the 
survey area JIdrs and within a radius r = R CO mp(My), 
we aim to determine a correction factor, c, that gives the 
total number of satellites brighter than My within a spher- 
ical volume associated with the larger outer radius -Routed 
Ntot = c(r, f2) iV bs(< r, < il). For our main results we 
count galaxies within -Router = 417 kpc, corresponding to 
the distance to the outermost Milky Way satellite (Leo T), al- 
though we do consider other i? utor values in Sj3] 

It is useful (although not precisely correct) to think of the 
correction factor as a multiplicative combination of a radial 
correction factor and an angular correction factor, c = c r CQ. 
The first correction, c r , will depend on the radial distribu- 
tion of satellites. If there is no systematic angular bias in the 
satellite distribution, and if the experiment is performed many 
times, then the second correction factor should have an aver- 
age value < cq >= 1//dr5 = 5.15. However, if the satel- 
lite distribution is anisotropic on the sky, the value of cq for 
any particular survey pointing may be significantly different 
from the average. An estimate of this anisotropy is essential 
in any attempt to provide a correction with meaningful errors. 
Below, we use the satellite halo distribution in Via Lactea to 
provide an estimate for the overall correction (see i j2.31 l. Note 
that in our final corrections presented in §3, we do not force 
the radial and angular corrections to be separable, but rather 
use a series of mock survey pointings within the simulation 
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FIG. 1. — The completeness radius for dwarf satellites. The three rising 
lines show the helio-centric distance, fJ conlp , out to which dwarf satellites 
of a given absolute magnitude are complete within the SDSS D R5 su rvey. 
The solid is for the published detection limits in Koposov et al. (2 007) and 
the other lines are the other two detection limits described in 812.21 The dot- 
ted black line at 417 kpc corresponds to our fiducial adopted outer edge of 
the Milky Way halo satellite population (-Router)- The data points are ob- 
served satellites of the Milky Way and the Local Group. The red circles are 
the SDSS-detected satellites, the only galaxies for which the detection lim- 
its actually apply, although the detection limits nevertheless also delineate 
the detection zone for more distant Local Group galaxies (yellow diamonds). 
Purple squares indicate classical Milky Way satellites. The faintest object 
(red triangle) is Segue 1 , which is outside the DR5 footprint. 
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2.2. SDSS Detection Limits 



Kop osov et al.l (l2007h constructed an automated pipeline to 
extract the locations of Milky Way satellites from the PR5 
stellar catalog. They then constructed artificial dSph galaxy 
stellar populations (assuming a Plummer distribution), added 
them to the catalog, and ran their pipeline. The detections 
of galaxies were use d to construct detec tion thresholds as a 
function of distance ( Ko posov et al.ll200~7l Figure 8). These 
thresholds are mostly constant as a function of surface bright- 
ness for /i < 30 and linear with the logarithm of distance. 
Hence, their detection threshold is well approximated as a log- 
linear relationship between the absolute magnitude of a dSph 
galaxy, My, and the volume, V comp (My), o ut to which the 
PR5 c ould detect it. Specifically, Figure 13 of Koposov et al. 
(120071) implies 
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such that the completeness volume follows Vc 0mp oc Ly 0//2 . 
We adopt this form as the dwarf galaxy detection limit of 
PR5. This volume may be related to a spherical completeness 
radius, -R C omp, beyond which a dSph of a particular magni- 
tude will go undetected: 

/ o \ 1/3 

( -4— 10(- aMv -^ 3 Mpc, (2) 

where /drs = 0.194 is the fraction of the sky covered by 
PR5. For our f i ducia l relation, we use the result presented in 
iKoposov et all d2007l Figure 13), which is fit by a = 0.60 
and b = 5.23. We note that for this value of a, we have 
#com P ~ 66 kpc(L/lOOOL ) 1/2 , and the SPSS is complete 
down to a fixed apparent luminosity. The implied relationship 
between galaxy luminosity and corresponding helio-centric 
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FIG. 2. — The cumulative count of Via Lactea subhalos as a function of 
the (current) maximum circular velocity of the subhalo (n ma x, solid), and 
as a function of the largest maximum circular velocity ever obtained by the 
subhalo (f pca ]j, dashed). 

completeness radius is shown by the solid line in FigureQ] We 
also consider two alternative pos s ibilitie s. One is obtained by 
fitting the data in Koposo vet al.l (12007) Table 3, which gives 
a = 0.684 and b = 5.667 (black dotted), and the other is 
a line that passes through the new SDSS satellites, on the 
assumption that some of them are of marginal detectability 
(a = 0.666 and b = 6.10; cyan dashed). 

It is important to recognize that, in principle, the detectabil- 
ity of satellites at a particular radius is not a step-function 
between detection and non-detection. It also should not be 
spherically symmetric (independent of latitude with respect to 
the disk plane), nor should it be independent of other variables 
(such as satellite color or background galaxy density). How- 
ever, iKoposov et al.l (120071) found that a simple radial depen- 
dence provided a good description of their simulation results, 
and we adopt it for our corrections here. They did, however, 
find that galaxies within this "completeness" boundary were 
not always detected at 100% efficiency (depending on their 
distance and luminosity) and published detection efficiencies 
for the known SDSS dwarfs (their Table 3). We use these 
published detection efficiencies in our fiducial corrections to 
the luminosity function below. We also investigate how our 
results change when we assume 100% efficiency in the cor- 
rection. 

For reference, the horizontal dotted line in Figure [TJ marks 
our fiducial adopted i? utcr radius for the Milky Way halo 
(slightly larger than the virial radius, in order to include Leo 
T). According to this estimate, only satellites brighter than 
My — —7 are observable out to this radius. The fact that 
the faintest dwarf satellite galaxies known are more than 4 
magnitudes fainter than this limit (see Table 1) immediately 
suggests that there are many more faint satellite galaxies yet 
to be discovered. 

2.3. Via Lactea 

iDiemand. Kuhlen. & Madaul (l2007al) describe the Via 
Lactea simulation, which is among the highest resolution 
ACDM N-body simulations of a Milky Way-like dark mat- 
ter halo yet published. The mass of Via Lactea is M200 — 
1.8 x 10 12 Mq, with a corresponding virial radius i?200 = 389 
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FIG. 3. — Inverse of cumulative subhalo counts within the specified radius. 
Here /(< r) is the fraction of subhalos that exist within a given radius, 
normalized to unity at -Router = 417 kpc. The radius is "heliocentric", 
defined relative to the (8,0,0) kpc position of the Via Lactea simulation. We 
include three populations of subhalos, as identified in the key. The inset 
focuses on the radial range where most of the new ultra-faint SDSS satellites 
have been detected. 

kpc, where M200 and -R200 are defined by the volume that 
contains 200 times the mean matter density. It resolves a large 
amount of substructure, recording 6553 subhalos with peak 
circular velocities larger than 5 km s _1 at some point in their 
history, of which 2686 are within our adopted i? utcr =417 
kp c at z = 0. We use the public data release kindly provided 
by IDie mand et al. (2007a) in what follows 3 . 

Figure 12 presents the cumulative maximum circular veloc- 
ity function, 7V~(> v max ), for Via Lactea subhalos with halo- 
centric radius R < 417 kpc at z = (solid), along with the 
cumulative "peak" circular velocity function (dashed) for the 
same halos 7Y(> t> p0 ak)- Here u pca k is the maximum cir- 
cular velocity tha t the subhalos ever ha d over their history. 
As emphasized bv lKravtsov et al.l (|2004) many subhalos have 
lost considerable mass over their history, and therefore w pea k 
may be a more reasonable variable to associate with satellite 
visibility than the (current) subhalo v max . The information 
contain ed within this figure is presented elsewhere in the lit- 
erature (Diema nd et al.l l2007a). although we include it here 
for the sake of completeness. 

An important ingredient in the luminosity bias correction is 
the assumed underlying radial distribution of satellites. We 
determine this distribution directly from Via Lactea. We are 
interested in determining the total number of satellites, N to t, 
given an observed number within a radius, N ^ s (< r). The 
correction will depend on the cumulative fraction of objects 
within a radius r compared to the total count within some 
fiducial outer radius: /(< r) = 7V(< r)/N tot . The associ- 
ated radial correction factor is then the inverse of the cumula- 
tive fraction, c r — / _1 (< r), such that 7V tot = c r AT b s (< r). 
This correction factor is shown in Figure [3] for three differ- 
ent choices of subhalo populations: u poa k > 10 km s _1 , 
Wmax > 7 km s _1 , and v pca k > 5 km s" 1 ( i.e. all subhalos 
in the Via Lactea catalog). 

In order to mimic a heliocentric radial distribution, in Fig- 
ure[3]we have placed the observer at a radius of 8 kpc within 

3 |http : / /www . ucolick . org/ ~{}diemand/vl/ data ■ html] 
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a fictional disk centered on Via Lactea. The orientation of 
the "disk" was chosen to be in the Via Lactea xy plane for 
all figures that show subhalo distributions, but we allow for 
a range of disk orientations for the full correction presented 
in $3] We do this for completeness, but stress that our re- 
sults are extremely insensitive to the choice of solar location, 
and are virtually identical if we simply adopt a vantage point 
from the center of Via Lactea. The total count is defined 
within our fiducial i? ou t i> such that /(< 417 kpc) = 1.0. 
An important result of this is that our correction does not 
depend on the number of subhalos in Via Lactea, only on 
the shape of the distribution. Still, the shap e varies among 
some sub-populations of subha lo. As noted inlKravtsoy et al.1 
(2004). lDiemand et al.l d2007al) and lMadau et al.l (120081) . sub- 
halos chosen to have a large u pca k tend to be more centrally 
concentrated. However, the correction we apply is fairly in- 
sensitive to the differences between these choices of subhalo 
populations. 

The most important corrections will be those for the faintest 
galaxies, which are just observable at local distances r < 50 
kpc. It is thus clear that the radial correction factors associ- 
ated with the three different subhalo populations show very 
little differences at the radii of relevance (see inset in Fig- 
ure [3). Motivated by this result, we use the full Via Lactea 
subhalo catalog in our fiducial model because it provides a 
larger statistical sample of subhalos. However, for the sake 
of completeness, we present the final corrected counts for our 
corrections using the w pca k > 10 km s _1 and v max > 7 km 
s _1 samples in Table |3]below. 

Via Lactea can also be used to extract the angular distri- 
bution of subhalos. This will be important for determining 
the errors on the overall satellite abundance that may arise 
from limited sky coverage. The SDSS DR5 footprint covers 
8000 square degrees, or a fraction Jdr5 = 0.194 of the sky. 
While this is a sizable fraction of the sky, the upper image in 
Figure |4] shows that the subhalo distribution (projected out to 
the virial radius, in this case) is quite anisotropic even on this 
scale. The color of each pixel is set by the fraction of all the 
subhalos, /(< ft), that exist within an area of fluR5 = 8000 
square degrees, centered on this pixel. We see clearly that 
the fraction of satellites within a DR5-size region can vary 
considerably (from 0.12 to 0.28) depending on the pointing 
orientation. The same information is shown as a histogram in 
the lower panel, now presented as the implied multiplicative 
"sky coverage" correction factor, co = l//(< Each of 
the 3096 pixels in the sky map is included as a single count in 
the histogram. Note that while the average correction factor 
is cn = 1 //oij5 = 5.15 (as expected) the correction can vary 
from just ~ 3.5 up to ~ 8.3 depending on the mock survey's 
orientation. 

Before moving on, we note that the expected anisotropy 
may contain a possible solution to the "missing inner satellites 
problem" described in Diemand et al. (2007a), which notes 
that there are ~ 20 subhalos in Via Lactea within the inner- 
most parts of the halo, while only a few galaxies are actually 
observed. While numerical effects may become important for 
these subhalos given their sizes and how close they are to the 
center of the host, it may also simply be an observational cov- 
erage effect. If only subhalos within 23 kpc (the distance to 
Segue 1) are considered, ~ 15% of the sky has one or zero 
subhalos visible in a survey the size of SDSS. This means that 
even if all of the subhalos host galaxies, there is a 15% chance 
that we will at most see a galaxy like Segue 1 . The fact that 
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FIG. 4. — Example of subhalo angular anisotropy. Upper: Hammer projec- 
tion map of the angular anisotropy in Via Lactea subhalos with i> pea k > 5 km 
s _1 . Colors at each point indicate the fraction of subhalos, /(< Q), within 
417 kpc that are contained within an angular cone of fioRS = 8000 square 
degrees, in order to match the area of SDSS DR5. Lower: Distribution of the 
angular correction factors, cq = l//(< ST), which must be applied to the 
count within 8000 square degrees in order to return the full number of sub- 
halos within 417 kpc. The histogram includes 1000 pointings, evenly spaced 
over the sky. For reference, the sky coverage of DR5 is }dR5 = 0.194, 
which implies an angular correction factor of X/foRS = 5.15. As it must, 
c = 5.15 matches the mean of the distribution, but not the median, which is 
c = 5.27. Note that these correction factors do not yet include the effects of 
radial incompleteness/luminosity bias. 

Segue 1 was only recently discovered shows that there is eas- 
ily room for a significant number of inner satellites that have 
gone as yet undetected simply because a survey like SDSS 
is necessary to uncover such objects, but over the entire sky 
instead of just a fifth. 

2.4. Observed Satellite Galaxies 

A variety of authors have identified dSph's in the DR5 foot- 
print, including some that straddle the traditional boundaries 
between globular clusters and dSph's. For example, Will- 
man 1 was originally of unclear classification (Willm ari et alJ 
l200lh . but is now gener ally recognized as being a dark 
matte r dominated dSph dMartin et all l2007t IStrigari et all 
2008). Table Q] lists the complete set of dwarf galaxies 
from DR5 used in this analys is. For the SPSS dwarfs, 
we use the luminosities from iMartin et al.l d2008l) . which 
presents a homogeneous a naly sis, for all of the SDSS dwarfs, 
as well as Ide Jong et alJ (120081) for Leo T, which is based 
on much deeper photometry. We note that these results 
are in generally good agreement with pre-existing analyses 
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TABLE 1 

Prop erties of known Milky Way Satellite gal a xies. Data are 
fromIbothun & Th ompson ( 1988); Mateo ( 1998); Gre bel et al.I 
d2003D : ISimon & GbhaH2007D : IMartin et al.I J2008T) ; Ide Jong et al.I 

J2008I) . 



Satellite M v L v [L e ] d SU n[kpc] i? ha lf[pc] a e b 



SDSS-discovered Satellites 



*Bo6tes I 


-6.3 


2.83 x 10 4 


60 


242 


1.0 


* Bootes II 


-2.7 


1.03 x 10 3 


43 


72 


0.2 


* Canes Venatici I 


-8.6 


2.36 x 10 5 


224 


565 


0.99 


* Canes Venatici II 


-4.9 


7.80 x 10 3 


151 


74 


0.47 


*Coma 


-4.1 


3.73 x 10 3 


44 


77 


0.97 


* Hercules 


-6.6 


3.73 x 10 4 


138 


330 


0.72 


*Leo IV 


-5.0 


8.55 x 10 3 


158 


116 


0.79 


*Leo T 


-8.0 


5.92 x 10 4 


417 


170 


0.76 


' Segue 1 


-1.5 


3.40 x 10 2 


23 


29 


1.0 


*Ursa Major I 


-5.5 


1.36 x 10 4 


106 


318 


0.56 


*Ursa Major II 


-4.2 


4.09 x 10 3 


32 


140 


0.78 


*Willman 1 


-2.7 


1.03 x 10 3 


38 


25 


0.99 



Classical (Pre-SDSS) Satellites 



Carina 


-9.4 


4.92 x 10 5 


94 


210 




* Draco 


-9.4 


4.92 x 10 5 


79 


180 


1.0 


Fornax 


-13.1 


1.49 x 10 7 


138 


460 




LMC 


-18.5 


2.15 x 10 9 


49 


2591 




Leo I 


-11.9 


4.92 x 10 6 


270 


215 


1.0 


* Leo II 


-10.1 


9.38 x 10 5 


205 


160 


1.0 


Ursa Minor 


-8.9 


1.49 x 10 5 


69 


200 




SMC 


-17.1 


5.92 x 10 s 


63 


1088 




Sculptor 


-9.8 


7.11 x 10 5 


88 


110 




Sextans 


-9.5 


5.40 x 10 5 


86 


335 




Sagittarius 


-15 


8.55 x 10 7 


28 


125 





a Satellite projected half light radius. b Detection efficiency from Kopos ov et aT] 
(2007).* Galaxy sits within the SDSS DR5 footprint. + Satellite is not used in 
fiducial LF correction. 



dWillman et all 120051: IBelokurov et al.|[2006t ISimon & Gehal 
120071: IWalsh et alJl2007allbh 7 Here, we also list the "classical" 
dwarfs, which were discovered before the SDSS, along with 
Segue 1, which was discovered from data for the SDSS-II 
SEGUE survey (IBelokurov et alJl2007l) . All of t he objects we 
list in this table have large mass-to-light ratios dMartin et alJ 
l2007t ISimon & Gehdl2007l:IStrigari et alJl2008h . 

For our fiducial co rrections, following the convention of 
Kopo sov" et alJ (120071) . we have not included Segue 1, as it 
does not lie inside the DR5 footprint and hence the published 
DR5 detection limits are not applicable. We do include Segue 
1 in an alternative correction scenario below (see Table 01. 
We do not correct the classical dwarf satellite galaxies for 
luminosity bias or sky coverage, because appropriate detec- 
tion limits for these classical dwarf satellites are unclear given 
that they are not part of a homogeneous survey like SDSS. 
We assume that all satellites within those magnitude bins 
would have been discovered anywhere in the sky, with the 
possible exception of objects at low Galactic latitudes, where 
Milky Way extinctio n and contamination become significant 
( Will man et al.l2004al) . This assumption is conservative in the 
sense that it will bias our total numerical estimate low, but it 
is only a minor effect, as our correction described in |J3] is 
dominated by low luminosity satellites. 

Before we use the radial distribution of Via Lactea subha- 
los to correct the observed luminosity function, it is impor- 
tant to investigate whether this assumption is even self consis- 
tent with the data we have on the radial distribution of known 



satellites. The relevant comparison is shown in Figure[5] We 
have normalized to an outer radius -Router = 417 kpc (slightly 
larger than the Via Lactea virial radius) in order to allow a 
comparison that includes the DR5 dwarf Leo T; this exten- 
sion is useful because the known dwarf satellite count is so 
low that even adding one satellite to the distribution increases 
the statistics significantly. 

The radial distribution of all 23 known Milky Way satel- 
lites is shown by the magenta dashed line in Figure [5] The 
four solid lines show radial distributions for four choices of 
subhalo populations: t he 65 largest y Dea k( upper) subhalos (65 
LBA) as discussed in iMadau etafl (120081) . u pea k > 10 km 
s _1 (upper-middle), t> poa k > 5 km s _1 (lower- middle), and 
« mal > 7 km s _1 (lower). We note that the all-observed pro- 
file is clearly more centrally concentrated than any of the the- 
oretical subhalo distributions. However, as shown in Figure 
[U our limited ability to detect faint satellite galaxies almost 
certainly biases the observed satellite population to be more 
centrally concentrated than the full population. 

If we include only the 11 satellites (excluding SMC and 
LMC) that are bright enough to be detected within 417 kpc 
(My < —7), we obtain the thick blue dashed line. This dis- 
tribution is significantly closer to all of the theoretical sub- 
halo distributions, and matches quite well within r ~ 50 kpc, 
where the incompleteness correction to the luminosity func- 
tion will matter most. It is still more centrally concentrated 
then the distribution of all subhalos, however, as has been 
noted in the past (at least for the classical satellites - e.g. 
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Willman et all l2004bt iDiemand et all l2004t iKravtsov et aT] 
2004 . In order to more rigorously determine whether the 
theoretical distribution is consistent with that of the 1 1 "com- 
plete" satellites, we have randomly determined the radial dis- 
tribution of fp ak > 5 km s _1 subhalos using 1000 sub- 
samples of 1 1 subhalos each. The error bars reflect 98 % 
c.l. ranges from these subsamples, although they are corre- 
lated and hence represent only a guide to the possible scatter 
about the points in the distribution. From this exercise, we 
can roughly conclude that the complete observed distribution 
is consistent with being a random subsample of the u poa k > 5 
km s _1 population. We also note that the -u poa k > 10 
km s _1 distribution fits even more closely, well within their 
68% distribution (errors on this thoretical distribution are not 
shown for figure clarity, but they are similar in magnitude to 
those on the central line). 

To further investigate the issue of consistency between the 
various radial distributions, we determine the KS-probability 
for various distributions. The results are shown in Table [2] 
where in each case the p-value gives the probability that the 
null hypothesis is correct (i.e. that the two distributions are 
drawn from the same underlying distribution). While the dis- 
tribution of all satellites is clearly inconsistent with the sub- 
halo distributions, choosing only the complete satellites im- 
proves the situation, giving a reasonable probability of the 
complete distribution matching the u pca k > 5 km s _1 dis- 
tribution, and an even better match for the distributions that 
are cut on v poa k, as expected from the error bars. 

Finally, we note that our definition of "complete" is con- 
servative, since no satellites fainter than My = —8.8 were 
known before the SDSS survey. If we include only the (five) 
satellites bright enough to be detectable to i? uter that ex- 
ist within the well understood area of the DR5 footprint, 
we obtain the radial profile shown by the dashed cyan line 
("DR5 Complete"), which is even less centrally concentrated 
than any of the theoretical lines. While this distribution is 
marginally consistent with the other distributions based on 
the KS-testp-values, caution is noted given the small number 
of data points (5). We conclude that while more data (from 
deeper and wider surveys) is absolutely necessary in order to 
determine the radial distribution of Milky Way satellites and 
to compare it with theoretical models, the assumption that the 
underlying distribution of satellites tracks that predicted for 
subhalos is currently consistent with the data. We adopt this 
assumption for our corrections below - specifically, we use 
the i> pea k > 5 km s _1 distribution for our fiducial case. While 
the i> pea k > 10 km s _1 distribution fits the observed satellites 
more closely, we show in Sj3] that the various subhalo distri- 
butions affect only the final satellite counts by a factor of at 
most ~ 1.4. We chose to use all Via Lactea subhalos (corre- 
sponding to w pca k > 5 km s ~ 1 ) as our fiducial case in order to 
reduce statistical noise from the smaller number of satellites 
in the w pea k > 10 population - note that numerical effects 
will generally produce an under-correction (giving a smaller 
number of satellites), as described in $3] 

Previous studies have pointed to discrepancies between 
simulations and observations of other galaxies' satel lite dis- 
tributions (e.g. lDiemand et al.ll2004tlChen et al.ll2006h . While 
this is of concern, much of these data are around clusters 
or other environments very different from that of the Local 
Group. Furthermore, the detection prospects for faint satel- 
lites around the Milky Way are of course far better than can 
be achieved for other galaxies, and the region of the lumi- 
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FIG. 5. — Radial distributions for various populations of observed satel- 
lites as compared to several sample subhalo distributions in Via Lactea. Solid 
lines are four populations of Via Lactea sub halos within 417 k pc: from top 
to bottom, 65 largest before accretion (see Madau et al. 2008, Figure 7), 
"peak > 10 km s" 1 , All (i.e. v pcak > 5 km s _1 ), and u max > 7 
km s _1 . The dashed lines are observed satellite distributions. From top to 
bottom, we have the"All Observed"(magenta), which consists of all known 
Milky Way dSph satellites; the "Complete" (blue) distribution corresponds 
to only those with magnitudes corresponding to fi comp > 417 kpc; "DR5 
Complete" (cyan) corresponds to satellites brighter than the completeness 
magnitude, and also in the DR5 footprint - these are the only observationally 
complete and homogeneous sample. The dotted line is the global dark matter 
distribution for the Via Lactea host halo. Error bars (98 % c.l. ) are de- 
rived by randomly sampling from the Via Lactea subhalos the same number 
of satellites as in the "Complete" sample (1 1). Note that this means the error 
bars only apply for comparing the solid lines to the dashed line, and that they 
are correlated, but still represent the scatter of individual bins in any possible 
observed sample of 1 1 satellites from Via Lactea. 



TABLE 2 
KS TEST RESULTS FOR RADIAL 
DISTRIBUTIONS SHOWN IN FIGURE^ 



Distribution 1 


Distribution 2 


p- value 


All Observed 


65 LBA 


57.6% 


All Observed 


"peak > 10 


10.3% 


All Observed 


All Subhalos 


0.4% 


All Observed 


^max ^> 7 


0.1% 


Complete 


65 LBA 


95.1% 


Complete 


"peak > 10 


38.2% 


Complete 


All Subhalos 


12.1% 


Complete 


^max ^ 7 


8.2% 


DR5 Complete 


65 LBA 


8.6% 


DR5 Complete 


"peak > 10 


14.4% 


DR5 Complete 


All Subhalos 


48.2% 


DR5 Complete 


"max ' 


63.1% 


All Observed 


Complete 


87.6% 


All Observed 


DR5 Complete 


6.6% 


Complete 


DR5 Complete 


41.1% 



nosity function that is of importance for this correction (e.g. 
My > —7) has never been studied outside the Local Group. 

3. LUMINOSITY FUNCTION CORRECTIONS 

The cumulative number of known Milky Way satellites 
brighter than a given magnitude is shown by the lower (red) 
dashed line with solid circles in Figure [6] This luminosity 
function includes both the classical dwarf satellites and the 
fainter, more recently-discovered SDSS satellites (excluding 
Segue 1 for reasons discussed above). By simply multiplying 
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the SDSS satellite count by the inverse of the DR5 sky frac- 
tion (1/ Jdrs) and adding this to the classical dwarf satellite 
count, we produce the (green) long-dashed line. This first- 
order correction provides an extremely conservative lower es- 
timate on the total Milky Way satellite count by ignoring the 
details of luminosity bias in the SDSS. 

We use the following series of mock surveys of the Via 
Lactea subhalo population in order to provide a more re- 
alistic correction. First, an observer is positioned at dis- 
tance of 8 kpc from the center of Via Lactea. We then de- 
fine an angular point on the sky from this location and use 
it to center a mock survey of solid angle Odrs = 8000 
square degrees. We allow this central survey position to 
vary over the full sky using 3096 pointings that evenly sam- 
ple the sky. Although we find that the absolute position of 
the observer does not affect our results significantly, we also 
allow the observer's position to vary over six specific lo- 
cations, at (x,y,z) = (±8,0,0),(0,±8,0), and (0,0, ±8) 
on the Via Lactea grid. We acknowledge that there are 
(contradictory) claims in the literature concerning whether 
satellite galaxies are preferentially oriented ( either parallel 
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20051 iMetzetal. 120071 


Kuhlen et alJl2007l IWang et alJl2008l 


Faltenbacher et al.l2008 


), and equally contradictory claims re- 



iBailin et al. 2005; Dutton et al. 2007)). If there were a prefer- 
ential orientation, then the appropriate sky-coverage correc- 
tion factors would need to be biased accordingly, but for our 
correction, we make no assumptions about the orientation of 
the "disk". Therefore, any uncertainty in the correct orienta- 
tion of the disk is contained within the errors we quote on our 
counts. In the end, we produce 18576 equally-likely mock 
surveys each with their own correction factors, and use these 
to correct the Milky Way satellite luminosity function for an- 
gular and radial incompleteness. 

For each of the mock surveys, we consider every DR5 satel- 
lite (i = 1, 11, sometimes 12 ) with a helio-centric distance 
within i?outcr = 417 kpc and determine the total number of 
objects of its luminosity that should be detectable. Specifi- 
cally, if satellite i has a magnitude My that is too faint to be 
detected at -Router ('•<?. if My > —7), then we determine the 
number of Via Lactea subhalos, N(r < i? comp , ft < OdrsX 
that are situated within an angular cone of size ^drs and 
within a helio-centric radius -R C omp(M). We then divide the 
total number of Via Lactea subhalos, iVtot by this "observed" 
count and obtain a corrected estimate for the total number 
satellites of magnitude M v : 

e = ^ o) 

N{r<R comp {M v ),n<fl DR5 )' 
If the SDSS satellite i is bright enough to be seen at iiouter, 
then -Rcomp is replaced by -Router in the above equation. In 
this case, the correction factor accounts only for angular in- 
completeness. This method allows us to produce distribu- 
tions of correction factors for each Afy of relevance. Note 
also that because the correction factor is the fraction of sub- 
halos in the cone, it depends only on the distributions of Via 
Lactea subhalos, rather than directly depending on the total 
number, rendering it insensitive to the overall subhalo count 
in Via Lactea. Furthermore, any subhalo distributions that 
do not have enough subhalos to fully accommodate all pos- 
sible pointings are simply given the correction factor l/7AT tot , 
meaning that numerical effects tend to unde rcorrect, produc- 
ing a conservative satellite estimate. 



Combined Correction 

• • Fixed Coverage Correction 

• •Uncorrected/Observed 




FIG. 6. — Luminosity function as observed (lower), corrected for only 
SDSS sky coverage (middle), and with all corrections included (upper). Note 
that the classical (pre-SDSS) satellites are uncorrected, while new satellites 
have the correction applied. The shaded error region corresponds to the 98% 
spread over our mock observation realizations. Segue 1 is not included in this 
correction. 

Three example distributions of number-count correction 
factors are shown in Figure [7] for hypothetical objects of lu- 
minosity My = —3,-5, —7 and corresponding completeness 
radii -R CO mp = 77, 194, and 486 kpc. We see that while the 
brightest objects typically have correction factors of order the 
inverse of the sky coverage fraction, ~ 5, the faintest objects 
can be under-counted by a factor of ~ 100 or more. Note that 
to produce helio-centric sky coverage maps, we must assume 
a plane in Via Lactea in which we consider the Sun to lie - for 
Figure [7] this is the xy plane, but for our final corrections we 
consider all possible orientations (described above). 

We construct corrected luminosity functions based on each 
of the 18576 mock surveys by generating a cumulative count 
of the observed satellites. We weigh each satellite i by its 
associated correction factor c 1 and (in our fiducial case) its 
detection efficiency e\ For each of the new satellites we use 
the quoted detection efficiencies from Kop osov et al.l (120071 
Table 3). We reproduce those efficiencies in our Table 1, and 
we assume e = 1 for all of the classical satellites that are not 
within the DR5 footprint. Explicitly, the cumulative luminos- 
ity function for a given pointing is: 



n{< M v ) 



<M V 



c(M v ) 



(4) 



For a given scenario (subhalo population, completeness lim- 
its, and detection efficiencies for each satellite) we determine 
the luminosity function for each pointing and disk orientation. 
We are then able to calculate a median luminosity function 
and scatter for each scenario. 

Our fiducial corrected luminosity function shown by the up- 
per blue solid line in Figure [6] and the shaded band spans 
the 49% tails of the distribution. Note that the errors van- 
ish around My = —9 because all satellites brighter than 
that are "classical" pre-SDSS satellites and are left com- 
pletely uncorrected on the conservative assumption that any 
objects brighter than this would have been detected previ- 
ously. Our fiducial scenario counts galaxies within a radius 
Pouter = 417 kpc and excludes Segue 1 from the list of cor- 
rected satellites because it is not within the DR5 footprint. 
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c(M v = -3) 




32 34 36 38 40 42 44 46 48 
LF correction factor 



c(M v =-5) 




9 10 11 12 13 
LF correction factor 



c(M v =-7) 




9.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 
LF correction factor 

FIG. 7. — Sky anisotropy maps and corresponding distributions for 3 sample 
limiting magnitudes. The sky projections are the corresponding correction 
factor (i.e. inverse of the enclosed fraction of total satellites) for a DR5-sized 
survey pointing in a given direction corresponding to a satellite of a particular 
magnitude. The histograms show the distribution of correction factors from 
evenly sampling the sky maps. From top to bottom, limiting magnitudes are 
M v = -3, -5, -7 corresponding to R CO mp = 77, 194, 486 kpc. 



In ad dition, this scenario u ses quoted detection efficiencies e 
from iKoposov et al.l (120071) and the completeness radius re- 
lation in Equation |2 with a = .6 and b = 5.23. With this 
fiducial scenario, we find that there are 398±j£ 8 (98% c.l.) 
satellites brighter than Boo II within 417 kpc of the Sun. 

We have performed the same exercise for a number of dif- 
ferent scenarios as described in Table [3] and summarized in 
Figure [8] Table [3] assigns each of these scenarios a number 
(Column 1) based on different assumptions that go into the 
correction. Scenario 1 is our fiducial case and counts all satel- 
lites brighter than Boo II {My > -2.7) within Router = 417 
kpc. Scenario 2 counts the total number of satellites brighter 
than Segue 1 (My > —1.5) if we consider Segue 1 as in- 
cluded in the DR5 sample (Column 5). Scenarios 3-4 exclude 
Boo II or Willman 1 (the faintest satellites) along with Segue 
1, considering the possibility that these satellites are low lu- 
minosity owing to strong environmental effects, or are not 
even truly dwarf galaxies at all. Alternatively, scenario 4 may 
be considered the total number brighter than Coma (the 4th 
faintest). Scenarios 5-10 allow for different choices within 
the Via Lactea subhalo population (Column 7), and 1 1-14 re- 
flect changes in outer Milky Way radius, i? utor (Column 4). 
Scenarios 16-17 allow different values for a and b in Equa- 
tion |2] (specified in Columns 2 and 3) for the R CO mp(Mv) 
relation, and 15 assumes that the detection efficiency for all 
satellites is e = 1. Finally, Scenario 18 makes the extreme 
assumption that there is no luminosity bias and includes only 
a sky-coverage correction factor. In columns 8-10, -/V sats is 
the number of satellites expected for a given scenario, while 
iV U ppe r and -/Vi owcl are the upper and lower limits for the 
98% distribution. The final column gives the faint-end (i.e. 
—2 > My > —7) slope if the luminosity function is ap- 
proximated by the form dN/dLy ~ L a . Note that some of 
these luminosity functions are not consistent with a power law 
within the anisotropy error bars, but we include the best fits 
for completeness (see $6]for further discussion). 

From Figure[S], for all of the cases that adopt the fiducial de- 
tection limits (scenarios 1-15), we may expect as many as 
~ 500 satellites within ~ 400 kpc. Even the most conser- 
vative completeness scenarios (4,10,11, and 14) suggest that 
~ 300 satellites may exist within the Milky Way's virial ra- 
dius. Scenario 17, which relies on a less conservative but not 
unreasonable detection limit, suggests that there may be more 
than ~ 1000 Milky Way satellites waiting to be discovered. 
We now turn to a discussion of the prospects for this exciting 
possibility. 

4. PROSPECTS FOR DISCOVERY 

Upcoming large-area sky surveys such as LSST , 
PES, PanSTARRS, and SkyMapper dlvezic et all 1200: 
The Dark Energy Survey Collaboration 120051 : iKaiser etal 



2002; Kell er et al.l l2007h will survey more sky and provide 
deeper maps of the Galactic environment than ever before 
possible. In this section we provide a first rough estimate for 
the number and type of satellite galaxies that one may expect 
to discover with these surveys. 

While it is impossible to truly ascertain the detection limits 
without detailed modelling and consideration of the sources 
of contamination for the magnitudes and colors these surveys 
will probe, the simplest approximation is to assume that all 
the characteristics of detectability are the same as those for 
SDSS aside from a deeper limiting magnitude. If we do this, 
we can estimate the corresponding completeness radius for 
each survey by the difference between the limiting 5 — a r- 
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TABLE 3 

The last four columns provide median, upper, and lower estimates (98 % range) for the Milky Way satellite count within 
router, as well as the faint-end slope (as the schechter a, i.e. dn/dl v oc l y ) for several different scenarios. see text for 

a description. 



Scenario 


a 


b 


■Router [kpc] 


Excluded Satellites 


Detection Eff. 


Subhalos 


N sats 


Nuppcr 


Nlowcr 


Faint-end slope(a) 


I 


-0.600 


-0 719 


417 


oegi 


Yes 


All 


398 


576 


304 


1 87 4- n 1 Q 
— i.o i zn u.-iy 


2 


-0 600 


-0 719 


417 




Yes 


All 


558 


88 1 


420 


_1 78 4- n 9fl 
— I.JO in u. 


3 


-0.600 


-0.719 


417 


Spirt RonIT 


Yes 


All 


328 


561 


235 


— 1.47 ± 0.18 


4 


-0.600 


-0.719 


417 


Sep1 RooII Wil l 


Yes 


All 


280 


505 


194 


— 1 fifi + 1 7 

j- . \jyj _1_ U.l 1 


5 


-0.600 


-0.719 


417 


Seel 


Yes 


finax ^ 5 


464 


749 


316 


— 1.92 ± 0.23 


6 


-0.600 


-0.719 


417 


Segl 


Yes 


fmax ^> 7 


427 


942 


258 


-1.90 ±0.32 


7 


-0.600 


-0.719 


417 


Segl 


Yes 


^peak ^ 


302 


508 


198 


-1.90 ±0.42 


8 


-0.600 


-0.719 


417 


Segl 


Yes 


Upeak > 10 


289 


521 


191 


-1.76 ±0.29 


9 


-0.600 


-0.719 


417 


Segl 


Yes 


"peak > 14 


260 


531 


167 


-1.72 ±0.32 


10 


-0.600 


-0.719 


417 


Segl 


Yes 


65 LBA 


224 


752 


128 


-1.67 ±0.44 


11 


-0.600 


-0.719 


200 


Segl 


Yes 


All 


229 


330 


176 


-1.70 ±0.18 


12 


-0.600 


-0.719 


300 


Segl 


Yes 


All 


322 


466 


246 


-1.80 ±0.19 


13 


-0.600 


-0.719 


389 


Segl 


Yes 


All 


382 


554 


292 


-1.87 ±0.19 


14 


-0.600 


-0.719 


500 


Segl 


Yes 


All 


439 


637 


335 


-1.89 ±0.20 


15 


-0.600 


-0.719 


417 


Segl 


No 


All 


184 


259 


142 


-1.65 ±0.19 


16 


-0.684 


-0.753 


417 


Segl 


Yes 


All 


509 


758 


381 


-1.95 ±0.20 


17 


-0.667 


-0.785 


417 


Segl 


Yes 


All 


1093 


2261 


746 


-2.15 ±0.26 


18 


N/A 


N/A 


417 


Segl 


N/A 


All 


69 


100 


50 


-1.16 ±0.21 
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Scenario 

FIG. 8 . — Number of Milky Way satellites for several different scenarios as 
described in Table|3] Error bars reflect 98% ranges about the median values 
(points). Scenario 1 (black circle) is our fiducial case and counts all satellites 
brighter than Boo II {My > -2.7) within .Router = AIT kpc. Scenario 2-4 
(blue triangles) consider inclusion or exclusion of different satellites. Sce- 
narios 5-10 (green squares) allow for different choices within the Via Lactea 
subhalo population, 11-14 (red diamonds) reflect changes in -Router, 15 (cyan 
hexagon) corresponds to no detection efficiency correction, and 16 and 17 
(magenta pentagons) are for different completeness limit assumptions. Sce- 
nario 18 (yellow cross) includes only the angular (sky coverage) correction 
factor, i.e. no luminosity bias correction. 

band point source magnitude for Sloan (Mgrjss = 22.2) and 
the corresponding limit for the new survey: 



plim 
^comp 

PSDSS 
l corap 



(5) 



Figure|9]shows the results of this exercise for several planned 
surveys. According to this estimate, LSST, with a co-added 
limit of A/lsst = 27.5, should be able to detect objects as 
faint as the faintest known galaxies out to the Milky Way virial 
radius. Moreover, LSST will be able to detect objects as faint 
as L ~ 1OOL (if they exist) out to distances of ~ 200 kpc. 

Examining the prospects for dwarf satellite detection with 
LSST in more detail, in Figure[l0]we plot the number of satel- 
lites that LSST would detect in Air sterradians assuming that 



the true satellite distribution follows that expected from our 
fiducial correction presented above. Remarkably, even with 
single-exposure coverage (to 24.5), LSST will be able to dis- 
cover a sizable fraction of the satellites in the Southern sky, 
> 100 galaxies by this estimate. Moreover, the co-added data 
would reveal all ~ 180 dwarf galaxies brighter than Boo II 
(and perhaps even more if they exist at lower luminosities). A 
survey such as this will be quite important for testing galaxy 
formation models, as even the shape of the satellite luminos- 
ity function is quite poorly constrained. In Figure [10] the 
Skymapper line appears close to a power law over the range 
shown, while the "actual" luminosity function is certainly not. 

Finally, in Table H] we present predictions for the number 
of satellites that will be found in each of a series of upcoming 
large-ar ea surveys, as w ell an example smaller-scale survey, 
RCS-2 (lYee et alj|2007l) . Given their limiting magnitudes and 
their fractional sky coverage and assuming the fiducial sce- 
nario outlined above, we can essentially perform the inverse 
of the correction described in Section [3] and determine the 
number of satellites the survey will detect (as done in Figure 
[TOlfor LSST). Note that this still ignores more complicated is- 
sues such as background galaxy contamination and problems 
with stellar crowding, which are potentially serious compli- 
cations, particularly for the faintest of satellites. Hence, these 
estimates are upper limits, but to first order they are a repre- 
sentative number of the possibly detectable satellites given the 
magnitude limits and sky coverage. Some of the footprints of 
these surveys will overlap, however, which will improve the 
odds of detection, but reduce the likely number of new detec- 
tions for later surveys (this effect is not included in the counts 
in Table©. 

5. DISCUSSION 
5.1. Caveats 

While our general expectation that there should be many 
more satellites seems to be quite solid, there are several im- 
portant caveats that limit our ability to firmly quote a num- 
ber or to discuss the expected overall shape of the luminosity 
function. 

1. Detection Limits. The simplifying approximation that 
Equation [T] describes the complexities of satellite de- 
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TABLE 4 

Predicted number of detectable satellites in a 
series of upcoming surveys (ignoring 
crowding/confusion effects, which will reduce the 
number of detections of faint obiects). 



Survey Name Area (deg 2 ) Limiting r [mag] N BS ,ts 



RCS-2 


1000 


24.8 


3-6 


DES 


5000 


24 


19-37 


Skymapper 


20000 


22.6 


42-79 


PanSTARRS 1 


30000 


22.7 


61-118 


LSST 1-exp 


20000 


24.5 


93-179 


LSST combined 


20000 


27.5 


145-283 



tection cannot be correct in detail - we expect a depen- 
dence on properties such as Galactic latitude and color. 
However, we are leaving these issues aside for this first- 
order correction. Mo reover, the detection limits from 
Kopo sov et al.l d2007l) apply only to objects with sur- 
face brightness fly < 30; if dSph's more diffuse than 
this exist, they will have been missed by SDSS. In this 
sense, our results are conservative, as even more very 
low surface brightness dwarf galaxies may be discov- 
ered by new, deeper surveys. 

2. Input Assumption. The assumption that the satellite 
population tracks the full underlying subhalo popula- 
tion is a major assumption of this correction. As Fig- 
ure[5]shows, this assumption is consistent with the cur- 
rent data, and modest cuts on v vea ^ improve the situa- 
tion further. Given the size of the error bars, the exis- 
tence of a discrepancy simply cannot be resolved with- 
out surveys that have better completeness limits. There 
are certainly reasons to expect that the satellite popu- 
lation will not track the subhalo population in an un- 
biased way - for example, if tidal forces play a role 
in making ultra-faint dSph galaxies. However, there 
are also reasons to suspect that there is no such bias. 
The dark matter subhalo masses of the known Milky 
Way dwarfs are approximately the same over a ~ 4 
orders-of-magnitud e spread in lumin osity (Strigari et 
al. 2008, in prep.; IStrigari et alll2008h . suggesting that 
the luminosity that any subhalo obtains is quite stochas- 
tic. Whether or not our input assumption about galax- 
ies tracing subhalos is correct, future surveys will pro- 
vide a means to test it, and thereby provide an important 
constraint on the formation processes of these extreme 
galaxies. 

3. Subhalo Distribution. Our correction relies on the Via 
Lactea subhalo population and its distribution with ra- 
dius and angle. While this provides a well-motivated 
starting point for this correction, Via Lactea is only a 
single realization of a particular mass halo and of a par- 
ticular set of cosmological parameters (WMAP 3 yr), 
and hence cannot necessarily be considered representa- 
tive of a typical halo in ACDM. Semi analytical models 
suggest that the cosmic variance in radial distributions 
should not be very large, as long as we consider ha- 
los that are small enough that the m ass-function is well 
populated (IZentner & Bullockll2003h . However, a suite 
of Via Lactea-type simulations will be necessary to pro- 
duce a distribution of corrections to statistically average 
over before this correction can be considered fully nu- 
merically robust. As discussed in Sjl] the uncertainty 
in the properties of the Milky Way halo also mean that 



using Via Lactea to accurately model the Milky Way 
is subject to those same uncertainties. This uncertainty 
is alleviated in this correction because only the fraction 
distribution matters - the absolute normalization is de- 
termined by the observed satellites rather than the sub- 
halo counts. But if it is true that the radial distribution 
of subhalos trac ks the global dark mat ter distribution, 
as suggested by Diema nd et alj d2007bl) and Figure [3] 
this uncertainty could be important if the inner portions 
of the radial distribution change substantially with cos- 
mological parameters or halo properties. 

Note that while the mass of Via Lactea is quite close 
to that obtained from ACDM-based mass models of th e 
Milky Way (e.g. lKlypin et alJ2002tlMadau et al.l2 008). 
the mass of the Milky Way halo is constrained by ob- 
servations only at the factor of ~ 2 level. The most 
relevant uncertainty associated with the possible mass 
difference between the Milky Way and Via Lactea is 
the difference associated with the uncertain radial scal- 

1/3 

ing. Since i?200 °^ M 2 q , a factor of ~ 2 mass uncer- 
tainty corresponds to a ~ ±25% uncertainty in virial 
radius. Fortunately, our corrections depend primarily 
on the relative radial distribution of satellites, not on 
their absolute masses, and the choice of exactly what 
radius within which galaxies are considered "satellites" 
is somewhat arbitrary, anyway. Moreover, the expected 
subhalo angular anisotropy (see fj3l> is more important 
for our corrections than the uncertainty associated in the 
Milky Way halo mass or precisely how the outer radius 
is defined. 

4. Input Luminosities/Satellites. Most of our corrections 
come from the faintest, nearby satellites, as they have 
very small R CO mp values compared to i? ou ter- Errors in 
the magnitudes of these faintest of satellites will result 
in significant changes in the derived luminosity func- 
tion. These errors are not propagated for this correction, 
but revisions to faint satellite magnitudes have the po- 
tential to significantly alter our estimates. Furthermore, 
the final count is highly sensitive to Poisson statistics in 
the innermost regions of the Milky Way, as well as any 
physical effect that biases the radial distribution of the 
faintest of the satellites (as discussed above). 

5.2. Implications 

Keeping in mind the caveats discussed above, it is worth 
considering the potential implications of a Milky Way halo 
filled with hundreds of satellite galaxies. For the sake of dis- 
cussion, we adopt A sat ~ 400, as obtained in our fiducial 
estimate. 

First, 7Y sat ~ 400 can be used to determine a potential char- 
acteristic velocity scale in the subhalo distribution. Under the 
(extreme) assumption that only halos larger than a given ve- 
locity host satellites, Figure [2] implies that N ~ 400 corre- 
sponds to a circular velocity v max > 7 km s _1 or a historical 
maximum of v pea k > 12 km s _1 . We note that these cutoffs 
correspond closely to those adopted in our corrections (Table 
0) and that there was no reason that they had to match. If we 
iteratively apply the correction in order to obtain perfect self- 
consistently between the corrected count and the input total, 
we find v pea k > 14 km s _1 . Interestingly, this scale is quite 
close to the thresh old where gas may be b oiled out of halos via 
photoevaporation dBarkana & Loebl[l999t) . We note however 
that the subhalo abundance at fixed mass or u max may vary 



Hundreds of Milky Way Satellites? 



11 



perha ps at as much as the fac tor of ~ 2 level from halo to 
halo (IZentner & Bull ock 2003). Therefore the "cosmic vari- 
ance error" will affect our ability to determine a characteristic 
subhalo mass based on counts (but does not affect the correc- 
tion itself, which depends only on the radial distribution, no t 
the total counts of subhalos). Recently, Die mand et al.l (120081) 
have published results from the Via Lactea II halo, which has 
a factor of ~ 1.7 times more subhalos at fixed v max than 
the Via Lactea halo we analyze here . If we take the veloc- 
ity function in Die mand et al.l d2008l) Figure 3, N sat would 
correspond to v max > 9 km s _1 . 

While the maximum circular velocity of a subhalo is a use- 
ful measure of its potential well depth, it is very difficult to 
measure v max direct ly from dwarf satellite stellar velocities 
(Strigari et al. 2007b). The best-determined dark halo observ- 
able is the integrate d mass within a fix ed radius within the 
stellar distribution dStrigari et al.ll200 7a). While the observed 
half-light radii vary from tens to hundreds of parsecs for the 
dwarf satellites, all dwarf satellites are found to have a com- 
mon mass scale of ~ 10 7 M Q within a fixed radius of 300 pc 
within their respective centers (Strigari et al 2008 in prep.) 
and to a similar extent a common mass of 10 6 within 100 pc 
(Strigari et al 2008). Although the masses within these scales 
are difficult to resolve with the Via Lactea simulation we con- 
sider in thi s paper, this mass sca le will be well-resolved in Via 
Lactea II (iDiemand et al.i 2008) and forthcoming simulations. 
In the future, a statistical sample of highly resolved subhalos 
will allow for a robust comparison between the dwarf satellite 
mass function and the subhalo population. This in turn will al- 
low corrections of the sort presented in this paper on various 
sub populations of subhalos and galaxies, providing a much 
more stringent consistency check between the mass function 
and luminosity function in ACDM-based models. It is im- 
portant to note, however, that the observed stellar kinemat- 
ics of the satellites does set stringent limits on their host halo 
Umax values. These limits are consistent with the results pre- 
sented here in that any of the reasonable sub-populations dis- 
cussed above (e.g. f max ^, 7 — 10 km s" 1 ) are n ot excluded 
by the current data (Stri gari et aTJ l2007bl 12008). We note 
that strong CD M priors suggest somewhat larger v max values 
(> 1 5km s _1 - IStrigari et ai1l2007bl [20081: iPefiarrubia et al l 
2008), which would (if anything) underpredict the observed 
visible satellite counts, according to our estimates. 

As we have shown, we expect the discovery of many 
more dwarfs to occur with planned surveys like LSST, DES, 
PanSTARRS, and SkyMapper. If so, it will provide important 
constraints on galaxy formation models, which, at present, are 
only poorly constrained by the present data. The ability to de- 
tect galaxies as faint as L ~ 100L Q provides an opportunity 
to discover if there is a low-luminosity threshold in galaxy for- 
mation, and to use these faint galaxies as laboratories to study 
galaxy formation in the extreme. The planned surveys will 
also provide an important measurement of the radial distribu- 
tion of faint satellites. This will help test our predictions, but 
more importantly provide constraints on more rigorous mod- 
els aimed at understanding why and how low-mass galaxies 
are so inefficient at converting their baryons into stars. 

Another important direction to consider is searches for 
similar satellites around M31. New satellites are rapidly 
being discovered in de ep surveys of its environs (e.g. 
McCo nnachie et all20 08). While there are indications of sub- 
stantial differences in the M3 1 satellites and their distribution s 
compared to the Milky Way (Mc Connachie & Irwinl 120061) . 
such comparisons are complicated by the fact that it is im- 




FlG. 9. — Maximum radius for detection of dSphs as a function of galaxy 
absolute magnitude for DR5 (assumed limiting r-band magnitude of 22.2) 
compared to a single exposure of LSST (24.5), co-added full LSST lifetime 
exposures (27.5), DES or one exposure from PanSTARRS (both 24), and the 
SkyMapper and associated Missing Satellites Survey (22.6). The data points 
are SDSS and classical satellites, as well as Local Group field galaxies. 



possible to detect the ultra-faint satellites that comprise most 
of our corrected satellite counts at the distance of M3 1 . with 
the much larger data samples that will be available with future 
deep surveys, it may be easier to compare the true luminosity 
function and distributions of Milky Way satellites to M3 1 and 
hence better understand the histories of both the galaxies, as 
well as better constrain how dSph's form in a wider ACDM 
context. 

Finally, if LSST and other surveys do discover the (full-sky) 
equivalent of ~ 400 or even ~ 1000 satellites, and appropriate 
kinematic follow-up with 30m-class telescopes like the Thirty 
Meter Telescope (TMT) confirms that these objects were in- 
deed dark-matter dominated, then it will provide a unique 
and powerful means to constrain the particle nature of dark 
matter. As discussed in the introduction, the mass function 
of dark matter subh alos is expected to rise steadily to smal l 
masses as ~ 1/M dKlvpin et al.lll999t|pTemand et alj|2007al) 
and the only scale that is expected to break this rise is the 
cutoff scale in the clustering of dark matter. When the MSP 
was originally formulated, scenarios like warm dark matter 
(WDM) were suggested as a means of "erasing" all but the 
~ 10 most massive subhalos per galaxy by truncating the 
power at ~ 10 8 M Q scales. If ~ 1000 satellites were dis- 
covered, then the same idea could be used to provide a limit 
on the small-scale clustering characteristics of the dark mat- 
ter particle. As a rough approximation, N ~ 1000 subha- 
los corresponds t o a minimum mass su bhalo in Via Lactea 
of M ~ 1O 7 M (IDiemand et al.ll2007al). or an original mas s 
(before infall) of M t ~ 3 x 10 7 M g (IDiemand et al.ll2007bl) . 
If we associate this z = subhalo mass with a limiting free- 
streaming mass, the n we obtain the bound m v J> 10 keV on 
the sterile neutrino ( Abazajian & Koushiappas 2006). This 
limit is competitive with the best constraints possible with the 
Lya forest and is not subject to the uncertainties of baryon 
physics. Of course, WDM simulations will be required in or- 
der to convincingly make a link between satellite counts and 
the small-scale power-spectrum, but these simulations are cer- 
tainly viable within the time frame of LSST. 
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FIG. 10. — Expected luminosity functions for LSST per Air sterradians. 
For the single exposure, the adopted limiting magnitude is r'i; m = 24.5, 
while for the co-added case, ri; m = 27.5 The Skymapper curve assumes 
a hypothetical full sky coverage survey of the same limiting magnitude as 
Skymapper (ri; m = 22.6). 

6. CONCLUSIONS 

The goal of this paper has been to provide reasonable, 
cosmologically-motivated corrections to the observed lumi- 
nosity function. Our primary aim is to motivate future 
searches for faint dwarf galaxies and to explore the status 
of the missing satellites problem in light of the most recent 
discoveries. By com bing completeness limits for the SDSS 
dKoposov et al.l2 007) with the spatial distribution of subhalos 
in Via Lactea, we have shown that there are likely between 
~ 300 and ~ 600 satellites brighter than Boo II within the 
Milky Way virial radius (Figure [5), and that the total count 
may be as large as ~ 2000, depending on assumptions (Table 
[3] Figure [8j. We also showed that the observed satellites are 
indeed consistent with tracing the radial distribution of subha- 
los, provided completeness limits are taken into account (Fig- 
ure [5] and Table |2). Moreover, we argued that future large 
sky surveys like LSST, DES, PanSTARRS, and SkyMapper 
should be able to see these satellites if they do exist, and 
thereby provide unprecedented constraints on the nature of 
galaxy formation in tiny halos. 

While this correction predicts that nearly all of the unde- 
tected satellite are faint (My > — 7) and consequently have 
low surface brightness, it is important to note that it is not 
clear how this maps onto satellites' w max . While a possible 
test of this correction's result lies in selecting sub-samples of 
the observed satellite pop ulation, f max is diffi cult to constrain 
in the known satellites dStrigari et alj|2007bh . and hence this 
exercise is suspect until simulations can resolve subhalos well 
enough to compare to observables such as the integrated mass 
within 300 pc (see jgj) . 

There are two major points to take away from this correc- 
tion: 

• As it was first formulated, the MSP referred to the mis- 
match between the then ~ 10 known dwarf satellite 

REFE] 



galaxies of the Milky Way and Andromeda, and the ex- 
pected couritofj^jrjo — 500 suWiak2s_with^ m a X > 10 
km s" 1 (!KlvpinetalJ[T999l:lMooreetal.lll999h . Our 
results suggest that the recent discoveries of ultra-faint 
dwarfs about the Milky Way are consistent with a to- 
tal population of ~ 500 satellites, once we take into 
account the completeness limits of the SDSS. In this 
sense, the primary worries associated with the MSP 
in CDM are alleviated. Nonetheless, it is critical that 
searches for these faint galaxies be undertaken, as the 
assumptions of this correction must be tested. 

• The shape of the faint end of the satellite luminosity 
function is not yet constrained well enough to deeply 
understand the theoretical implications. There exists 
yet a large parameter space in galaxy formation theory 
that will fall inside the error bars of Figure |6j and an 
even larger parameter space of models that are viable 
if our caveats and scenario possibilities are considered. 
Our results strongly suggest that the luminosity func- 
tion continues to rise to the faint end, with a faint end 
slope in our fiducial scenario given by dN/dMy — 

10 (0.35±0.08)M v +3.42±0.35 5 Qr dN / dLy K L « with 

a = — 1. 9 ± 0.2. This i s substantially steeper than the 
result of iKoposov et alj d2007l) . possibly the result of 
using a ACDM-motivated subhalo distribution instead 
of the analytic profiles used in that paper. But caution 
is advised in reading anything into the details of the 
shape, as nearly anything could be hiding within the 
faintest few bins when all the scenarios are considered. 
This is apparent from Figure[l0l where the Skymapper 
line appears as a power law over the range shown, while 
the actual luminosity function is certainly not. Further- 
more, depending on which scenario is used, the slope 
(a) can vary anywhere from -1.16 to -2.15 (see last col- 
umn of Table O. 

Fortunately, future deep large sky surveys will detect very 
faint satellites out to much larger distances and hence firmly 
observe the complete luminosity function out beyond the 
Milky Way virial radius (see Section [4). With these data, it 
will be possible to provide stringent limits both on cosmology 
and galaxy formation scenarios (see Section [5T2T >. Nonethe- 
less, the current data are not deep enough, and until the new 
survey data are available, there will be no way to put the spec- 
tre of the MSP completely to rest. 
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